Higher-order Kerr effect and harmonic cascading in gases 
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The higher-order Kerr effect (HOKE) has been recently advocated to explain measurements of the 
saturation of the nonlinear refractive index in gases. Here we show that cascaded third-harmonic generation 
results in an effective fifth order nonlincarity that is negative and significant. Higher-order harmonic cascading 
will also occur from the HOKE, and the cascading contributions may significantly modify the observed 
nonlinear index change. At lower wavelengths cascading increases the HOKE saturation intensity, while for 
longer wavelengths cascading will decrease the HOKE saturation intensity. © 2012 Optical Society of America 
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In homogeneous, isotropic media the nonlinear vari- 
ation of the refractive index can be expressed in terms 
of a power series of the light intensity /: An = nil + 



n 4 P 



-n^m.!™ 1 - However, higher-order nonlinearities 



have usually been regarded as negligible, as compared 
to the leading third-order nonlinearity ni (Kerr effect). 
A recent experiment aimed at characterizing the higher- 
order Kerr effect (HOKE) [Q has sparked a vivid debate 
on the actual role of these higher-order terms in non- 
linear beam-propagation dynamics. Indeed, the effective 
saturation of the HOKE nonlincarity could lead to fil- 
amcntation of femtosecond pulses in gases without the 
plasma playing an active role [0. The mechanism advo- 
cated in this case to arrest collapse of the wave packet 
was the saturation of the nonlinearity rather than the 
formation of plasma, as commonly accepted ||. Exper- 
iments have either confirmed [[|-((| or ruled out f J a 
preponderant role of HOKE in filamentation dynamics. 

In the strong cascading limit cascading of harmonics 
gives effective higher-order nonlinearities that compete 
with the inherent material ones. 2nd-order \^ '■ X^ cas- 
cading gives an effective x^ nonlinearity |l(],[ll[] , which 
can be made ultrafast and negative (self-defocusing) and 
even strong enough to overcome the self-focusing mate- 
rial Kerr nonlinearity [ jT2| . Similarly, 3rd-ordcr x^ '■ X^ 
cascading gives an effective x^ (0J- Here we study 
cascading contributions to the nonlinear HOKE coef- 
ficients in gases. This is unlike previous investigations 
that focussed on harmonic yield from cascading [||-^). 
Cascading relies on a phase mismatch Aki m = k(uj m ) — 
mk(w\) ^ between the fundamental wave (FW, u>±) 
and the mth harmonic (w m ). After two coherence lengths 
27r/|Afci m | back-conversion to the FW is complete, and 
as Afci m ^ the harmonic has a different phase veloc- 
ity than the FW. Thus, the back-converted FW photons 
are phase shifted with respect to the unconverted FW 
photons. For a strong phase mismatch Ak\ m L ^> 1 (L is 
the interaction length) , this process repeats many times 
and the FW feels a nonlinear phase shift parameterized 
by an effective higher-order nonlinear index n2 m , C asc oc 



_^(™)] 2 /A/ci m . We show that in gases cascading can af- 
fect the nonlinear index change and in particular the sat- 
uration intensity. Specifically, to lowest order cascading 
comes from a phase-mismatched third-harmonic genera- 
tion (THG) process, which we show generates an ri4 term 
that is negative and of significant strength. Similarly, 
cascading of HOKE terms gives contributions to sub- 
sequent orders of nonlinearity (rig, ng,,...) with values 
comparable with recent experimental Q and theoretical 
HOKE coefficients This implies that in the strong 
cascading limit, defined later, cascading may contribute 
significantly to the nonlinear dynamics. 

We now discuss cascaded THG, and will generalize to 
higher harmonics later. In the slowly-varying envelope 
approximation the inks electric-field coupled-wave prop- 
agation equations for THG u>\ + ui\ + ui\ — > lo 3 in the 
plane- wave limit in an isotropic lossless medium are 

(id z - ifc 2 (o;i)0 rT ) Ei 

' + 2E 1 \E 3 \ 2 + (El) 2 E 3 e lAkl3Z ] = (1) 

(id z - idi 3 d T - \k 2 {u 3 )d TT ) E 3 

+T& [E 3 \E 3 \ 2 + 2E 3 \E 1 \ 2 + \Eler^ z ] = (2) 

fc(w) = n{uj)uj / c is the wave number, n(u) the linear 
refractive index, k m (uj) = d m k(uj) / duj m the dispersion 
coefficients, d\ 3 = fci(oji) — fci(w3) the group-velocity 

mismatch, and T^) = 3u)jX^ /Scn(uij), j = 1,3. 
In the strong cascading limit, Afc 13 L ;§> 1 and no 
FW depletion are assumed Jl7|], so third harmonic 
(TH) self- and cross-phase modulation in Eq. (||) arc 
negligible and only the E^/3 term is relevant. The 
ansatz E 3tCa _ sc (z, r) = e _lAfel3Z -03( T ) makes Eq. (||) a 
temporal ODE. Neglecting dispersion for the moment 



[E^l 



gives E 3casc (z,T) 



-iAk 



>*E*(z,T)T™/(3Ak 13 ): 



the TH is effectively slaved to the FW. Inserting 
this result in Eq. (hi) yields the nonlinear term 



ifMl^il 2 \1 - l^r&VAftxa + 2|£ 1 | 4 (ir i 7Afc 13 ) 

where we used r^/3 ~ an d n(u}i)/n(uj 3 ) ~ 

Converting to intensity E\ — > [2/ ' Eon^ijc] 1 / 2 A\ 



I2r( 3 ) 



(3) 



1 



Tabic 1. Material M and cascading nonlinearitics at Ai = 0.8 /xm, p = 1 atm and T = 293 K. 





Afc 13 a 


ri2 


n4 


7^4, case 


n 6 


^6, case 


n 8 


^8, case 


n-10 ™10,casc 


Ar 


442 


1.0E-23 


-3.6E-42 


-1.8E-42 


40E-59 


0.2E-59 


-17.2E-76 


-3.0E-76 


0.9E-93 2.0E-93 


Air 


482 


1.2E-23 


-15.0E-42 


-2.3E-42 


21E-59 


0.8E-59 


-8.0E-76 


-1.8E-76 


1.0E-93 


o 2 


772 


1.6E-23 


-51.6E-42 


-2.6E-42 


48E-59 


2.2E-59 


-21.0E-76 


-3.3E-76 


2.2E-93 


N 2 


501 


1.1E-23 


-5.3E-42 


-1.9E-42 


14E-59 


0.3E-59 


-4.4E-76 


-1.0E-76 


0.5E-93 



"Calculated using Sellmeier equations |15||l6|. All units are SI (i.e. [A/C13] = m 1 . [ri2] = m 2 /W etc.). Note: the cascading calculations 



gave orders beyond nio, C asc, which are included in Fig. |lj. The propagated error in the leading cascading term ri4 iCasc is 12-35%. 



get (naAil^il 2 + n 4 ,ca SC Ai|Ai| 4 + r^LcAl^ilVi/c, 
where 77,2 = 3x /^ £ octi 2 (uji) is the Kerr nonlinearity. 
Thus, the cascaded TH gives rise to 5th- and 7th-order 
self-action nonlinearitics, where the 5th-order term is 

ri2 u>i 



TI4,, 



Ah 



3 C 



(3) 



This adds to the material quintic nonlinearity 77,4 so 
'Htot = n 4 + "4.casc- In other words, when THG is 
included in the description, '■ X cascading leads 
to a 5th-order HOKE term that is tunable in sign and 
magnitude. In the case we consider, a nonlinear medium 
without waveguide contributions pumped far from reso- 
nances, Afci3 > 0, so ri4. casc is always negative. 



At the next order we get 7i6,, 



(3) 

la 
o.casc 



+ n 



(5) 

6, case 



The former term was derived above, and the latter term 
comes from cascaded THG in the quintic nonlinear FW 
term (governed by 714), i.e. x'' 3 '' : X cascading. These 
add to the material n§, so tot = uq + uq casc , and read 



o,casc 



4 5) 

o.casc 



5rt2^4 i^i 
2Ak 



13 C 



(4) 



The expressions for the following orders are lengthy, 
but briefly in the strong cascading limit the harmonic is 
reduced to E m ^ casc oc — (Akim)' 1 E™e~ lAklmZ . Inserting 
these cascading terms in the FW equation, which con- 
tains all possible harmonic self- and cross-phase modu- 
lation terms, gives the cascaded higher-order self-action 
nonlinear coefficients. In this way we calculated HOKE 
contributions from 3rd- to llth-harmonic generation 
(nonlinearitics governed by 712 to nio), and Table 1 shows 
the cascading terms up to nio jC asc for various gases at 
Ai = 0.8 /urn. Note firstly that the cascading terms have 
similar strengths as the material nonlinearities at that 
level. Secondly, the Afci3 values for THG imply that the 
strong cascading limit can be achieved with propagation 
lengths fractions of a meter (the higher-order Aki, n are 
always larger than A/ci3, so they are less critical). 

Figure [j](a) shows the impact of cascading at Ai = 
0.8 yum, where we plot the variation of the nonlinear re- 
fractive index as a function of the pump intensity for 
the gases in Table [l| It is clear that the inclusion of the 
cascading effect shifts the inversion intensity 7i nv , where 
An = 0, to higher values. However, at longer pump wave- 
lengths the flattened gas dispersion enhances cascading. 
As a result the inversion intensity becomes smaller and 
departures from the ideal Kerr effect is observable al- 
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Fig. 1. (Color online) (a) The nonlinear index change Art 
vs. pump intensity for Ai = 0.8 /im, (b) the inversion 
intensity 7i nv vs. FW wavelength. In both plots the thin 
lines are the data from 0, while the bold lines include 
cascading corrections. In (b) the nonlinear coefficients 
include frequency dispersion using Miller's formula [hsl. 



ready at low intensities. In Fig. [l](b) the inversion intensi- 
ties are plotted vs. pump wavelength to summarize these 
trends. Without cascading, only a slight wavelength scal- 
ing [^8| of the nonlinearities is seen. When including cas- 
cading the inversion intensity peaks in the near-IR, and 
is pushed to become lower than without cascading at the 
beginning of the mid-IR. The calculations are based on 
the HOKE material nonlinearities from |pj, b ut we found 
similar results using the coefficients from [p.4[ (seemingly, 
the HOKE coefficients of [0 and [Q do not match, but 
this is in part due to the large difference in the num- 
ber of HOKE orders). Note that the diatomic gases (air, 
O2 and N2) also have fractional contributions from rota- 
tional and vibrational Raman effects that contribute to 
the overall observed nonlinearity as well, and this might 



2 



E%r L ] > 1 the 
,L ~ 1 



make it more difficult to observe effects of cascading. 

The strong cascading limit [jlTj assumes that the co- 
herence length Z co h = 7r/|Afc lm | is much smaller than (a) 
the interaction length (Afc lm L ^> 1) and (b) the char- 
acteristic nonlinear length (Afc im / '[ri™' E^r 1 ] > 1). To 
lowest order no pump depiction is found: i?3 iCasc above 
Eq. (||) is just phase modulated in z. When Aki m L ;§> 1 
cascading gives an effective Kerr-like nonlinearity (An 
independent on z), and when Afc lm /[r^ 
approximation An = n<2,jl 3 holds. For Ak\ 
this breaks down as the cascaded harmonic becomes am- 
plitude modulated by zsinc(Afci rn z/2). The cascading 
terms must then be corrected by 1— sinc(Afci m L) (resolv- 
ing the unphysical divergences at Ak\ m = 0). Nonethe- 
less, even for Aki m L < 1, strong cascading is possible if 
Akim/iT^E^r 1 } < 1 (strong FW depletion), but the 
cascading will then grow step- wise in z [jl) , so the cas- 
cading nonlinear coefficients cannot be parameterized. 

The increased cascading contributions predicted at 
longer wavelengths in Fig. ^(b) rely on reduced phase- 
mismatch values. Thus, fulfilling Aki^L ^ 1 can be- 
come problematic except for very long cells. Filamenta- 
tion experiments will therefore make an ideal test bed 
for cascading as they inherently use very long interac- 
tion lengths, where the sinc(Afci3L)-oscillations vanish. 
Alternatively, a solution could be to increase the pres- 
sure: the inversion intensity remains the same, but the 
phase mismatch values increase and overall we find that 
the cascading contributions to An are enhanced. 

The experiment behind the HOKE coefficients 0] used 
a pump-probe interaction length L 2± 1 mm |l9|| . Since 
the leading phase-mismatch values (see Table |I|) are 
around 1/mm or less, the cascading contribution is re- 
duced because the 1 — sinc(Afci3L) term must be applied. 
Specifically, we found that cascading corrections to all 
HOKE terms for the various gases are within the error 
bars of the measurements. Thus, the HOKE values in JlJ 
correspond to the intrinsic material nonlinearities. 

Even if cascading relies on a chain of up- and down- 
conversion steps, it can be considered instantaneous, 
much like an electronic Kerr nonlinearity: with full dis- 
persion, the cascading term in frequency domain be- 
comes E 3tCasc (z,il) oc -Afc 13 1 J r [£'j i (z,T)]i? casc (^) where 
J-[-] denotes the Fourier transform and R casc (fl) = 
Aki3/[k(il+ujs) — ki(u!i)il~3k(u>i)] is the cascading "re- 
sponse" function |l7j] . Neglecting higher-order dispersion 
the criterion df 3 < 2fc2(w3)Afci3 approximately deter- 
mines when the cascading response is nonresonant and 
thus broadband |T^]. We checked that this inequality 
is always fulfilled, and moreover since octave-spanning 
band widths are supported the corresponding temporal 
response of the cascading is extremely fast [ |12[ . 

Concluding, cascading in Kerr media from phase- 
mismatched harmonic generation may significantly con- 
tribute to the observed nonlinear index change in gases. 
The calculations rely on the strong phase-mismatch 
limit, and the lowest order contribution from 712, where 



^,(3) . ^,(3) cascading from a phase-mismatched third- 
harmonic gives an effective term, is always negative 
and significant compared to the material n^. A unidirec- 
tional pulse propagation equation can model such har- 
monic dynamics with a proper nonlinear expansion [ ^o| . 
For higher-order harmonics the conclusions were similar; 
such terms rely on including higher-order material non- 
linearities. At 800 nm the cascading from ??2 and higher- 
order coefficient (n^ up to nio) delayed the saturation 
point to higher intensities than what the material non- 
linearities dictate. This might explain why indications 
of saturation for near-IR filaments have been elusive. In- 
stead for longer (mid-IR) wavelengths, cascading pushed 
the saturation point to lower intensities, which will en- 
hance saturation and possibly affect filamcntation. 
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